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The properties of the superconducting and the anomalous normal state have been described 
by using the Eliashberg method. The pairing mechanism has been reproduced with help of the 
Hamiltonian, which models the electron-phonon and electron-electron-phonon interaction (EEPh). 
The set of the Eliashberg equations, which determines the order parameter function (tp), the wave 
function renormalization factor (Z), and the energy shift function (x) has been derived. It has 
been proven that for the sufficiently large values of EEPh potential, the doping dependence of order 
parameter (ifilZ) has the analogous course to that observed experimentally in cuprates. The energy 
gap in the electron density of states is induced by Z and x ■ Ike contribution from ip is negligible. 
The electron density of states possesses the characteristic asymmetric form and the pseudogap is 
observed above the critical temperature. 
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I. INTRODUCTION 

Superconductors, in which the condensate of the 
Cooper pairs is induced by the electron-phonon (EPh) 
interaction, can be modeled by Frohlich Hamiltonian dl¬ 
ls]. Starting from Frohlich operator, and using the canon¬ 
ical transformation to eliminate the phonon degrees of 
freedom, the effective Hamiltonian of the BCS theory 
can be obtained n mi- Then, the equation for the 
order parameter is derived. It should be noted that, 
due to the used approximations, the present model pro¬ 
vides the quantitative description of the superconducting 
state only in the weak-coupling limit, so in the fairly nar¬ 
row group of the low-temperature superconductors lanj. 
An important alternative to the above scheme is the ap¬ 
proach taken by Eliashberg [3]. In the case under con¬ 
sideration, the system of equations determining the ther¬ 
modynamic properties of the superconducting state is ob¬ 
tained directly from the Frohlich Hamiltonian. Note that 
the required calculations are performed using the Green 
matrix functions. Then, the Dyson formula is obtained at 
the level of the second order equations of motion isnn]. 

The Eliashberg equations can be divided into three 
groups. The first one allows to determine the function 
of the order parameter ( 99 k (*w„)), the second one serves 
to calculate the value of the wave function renormaliza¬ 
tion factor (Zk (iojn)), and the third one determines the 
energy shift function (xk (*w„)). The Eliashberg system 
is supplemented by the equation for the chemical poten¬ 
tial. The quantities k and denote respectively the mo- 
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mentum of the electron and the fermion Matsubara fre¬ 
quency: ojn = (tt//!) (2n — 1), where f3 = {kBT)~^ (fcs is 
the Boltzmann constant). It should be noted that the or¬ 
der parameter in the Eliashberg formalism is defined by: 
Ak(*w„) = (p]^{iLUn) The function Zk 

determines the ratio of the electron effective mass to the 
electron band mass, and the quantity Xk (*Wn) directly 
renormalizes the electron band energy. The great advan¬ 
tage of the approach proposed by Eliashberg is the fact 
that the thermodynamic properties of the superconduct¬ 
ing state can be determined under this scheme for any 
value of the electron-phonon coupling constant. The re¬ 
sults obtained for the classical superconductors usually 
provide the quantitative agreement between the theoret¬ 
ical predictions and the experimental data [zunHH]. 

The Eliashberg equations can be reduced to the BCS 
model, but it is necessary to skip the retardation, strong¬ 
coupling, and many-body effects: Ak (*Wn) Ak, 
Zk = 1, and Xk (*w„) = 0. 

The main purpose of the presented paper is to analyze 
the properties of the energy gap in the electron density 
of states in the framework of the Eliashberg formalism. 
We are interested in the case of the gap, which exists also 
above the critical temperature (Tc). Such gap is usually 
described as the pseudogap in difference of the energy 
gap induced by the order parameter. The inspiration for 
the conducted research was the experimental data ob¬ 
tained in the high-temperature superconductors based on 
cooper HMD- In the considered group of the materials, 
the pseudogap is observed, inter alia, in the photoemis¬ 
sion spectra [18], by the tunneling spectroscopy [TM2T] . 
and in the ARPES data However, we do not 

set ourselves to quantitative characterize the pseudogap 
in cuprates due to the fact of neqlecting of the Coulomb 
electron correlations, which are playing the important 
role in the considered superconductors pSIt28j . 
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II. ELIASHBERG APPROACH, PSEUDOGAP, 
AND PAIRING MECHANISM 


We demand that the pseudogap would appear sponta¬ 
neously in the electron density of states, and not be arti¬ 
ficially postulated in the formalism. In order to initially 
approximate the considered issue let us take into account 
the isotropic spectral function (A(a;,e)) that uniquely 
determines the electron density of states: N (w) = 
A(a;,0). The quantity A(w,e) should be calculated on 
the basis of the formula: A {uj,e) = —^Im [Gdiag (w,£)], 
while the diagonal part of the Green function takes the 
form: 


Gdiag : £) 


Z {u;)oj + e + xi^) 

{Z (w) uf -{e + x (w))^ - (w) 


( 1 ) 


Let us note that on the real axis Z (oj), x (w), and (p (w) 
take the complex values, thus the formula Q can de¬ 
note a very complicated function. For the case T > Tc, 
we assume p (w) = 0. However, the disappearance of 
the physical value of the order parameter above the criti¬ 
cal temperature (A (T) = p (T) jZ (T)), determined with 
the help of the equation: 


A(T)=Re[A(a; = A(T),r)], (2) 


does not automatically mean the disappearance of the 
function p (oj) on the whole real axis. Mindful of this no¬ 
tice, we write down the simplified formula on the electron 
density of states: 


N (oj) 


Zj (uj) uj-xi (w) 


(Zn (oj)oj-e-xi ( w ))^ -f (Zj (oj) oj - xi ( w ))^ 

(3) 

The indexes R and / denote the real and imaginary part, 
respectively. On the basis of the equation ([^, we con¬ 
clude that if the pseudogap in the electron density of 
states has to be observed above Tc, is will has to be 
induced by the strong-coupling and many-body effects 
(Z ( oj ) and x(^))- Hence, the following conclusions can 
be drawn: 

- The pseudogap cannot be described in the mean-field 
approximation, where Z (oj) = 1 and x (oj) = 0. It is 
the dynamic effect, i.e. it is related to the dependence of 
the wave function renormalization factor and the energy 
shift function on the frequency. The description of the 
pseudogap in the static approach is incorrect. 

- The pseudogap has nothing in common with the ad¬ 
ditional electron phase, because Z (oj) and x do not 
represent the order parameter. 

- Z (oj) and x (w) take the non-zero values both below 
and above Tc- Thus, they may make the important con¬ 
tribution in the energy gap in both superconducting and 
’’normal” state. 

- Due to the weak dependence of Z (oj) and x (‘j-’) on 
the temperature, in the case of the significant pseudo¬ 
gaps existence, we can expect the anomalous tempera¬ 
ture dependence of the energy gap in the superconduct¬ 
ing state. The energy gap in the superconducting state 


should smoothly transit into the pseudogap. As it was al¬ 
ready mentioned, our goal is to quite simply describe the 
properties of the energy gap existing also above Tq- For 
this reason, we consider incomplete pairing mechanism 
associated with the superconducting state in cuprates. 
We limit to the electron-phonon component, which con¬ 
sists of the electron-phonon and electron-electron-phonon 
ineraction. This choice was based on the following exper¬ 
imental and theoretical facts: 

- (i) in cuprates, the ARPES method clearly indicates 
the existence of the low-energy kink in the energy spec¬ 
trum near the phonon energy [551 US] , (h) with its help 
one can also observe the isotopic effect associated with 
the real part of the self-energy [30| . Although, these data 
are strongly undermined in the paper m, (iii) the iso¬ 
tope effect is seen for the critical temperature, especially 
in the highly undoped area |32|, (iv) the direct obser¬ 
vations of phonons can be obtained by using the scan¬ 
ning tunneling microscopy [53], and (v) phonons affect 
the penetration depth and Raman measurements [34] . 
[35] . From the theoretical point of view, the idea to use 
the Hamiltonian that explicitly takes into account the 
phonon degrees of freedom for the description of the high- 
temperature superconducting state is not new. The first 
advanced attempts in this direction were performed by 
Kim and Tesanovic [35] . The authors shown that strong 
Coulomb correlations may not completely suppress the 
phonon mechanism in the large doping regime. 

- The important role of EEPh interaction in cuprates 

is evidenced by the experimentally observed half-value 
of the magnetic flux It suggests that the electron 

quartets are formed. The existence of quartets can be the 
most easily reproduced with the help of the EEPh inter¬ 
action - for this purpose is sufficient to use the canoni¬ 
cal transformation that eliminates the phonon degrees of 
freedom [381 El] ■ It is worth to emphasize that ability to 
complete the BCS theory with the four-fermion interac¬ 
tion was first time noted by Rickayzen [40] . who moti¬ 
vated his idea as analogous to nuclear physics, where al¬ 
pha particle represents the stable system of four fermions. 
In the context of the description of the superconducting 
state, Rickayzen’s idea was followed after many years in 

[HHis]. 

- Carried out so far theoretical studies suggest that us¬ 
ing the Hamiltonian including EPh and EEPh terms it 
is possible to reproduce some properties of the supercon¬ 
ducting state in cuprates. In particular, it was shown 
that the model allows to correctly calculate the critical 
temperature and the energy gap at zero Kelvin. It should 
be noted that appropriate calculations were carried out 
for the large number of the chemical compounds, and 
each time obtained results were consistent with the ex¬ 
perimental data [33] [33] I33H33] . 

HI. THE FORMALISM 


The Hamiltonian, which models the electron-phonon 
and electron-electron-phonon interaction, in momentum 
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representation has the following form [351ISS] : 

(4) 

where first term describes the system of non-interacting 
electrons and phonons: 

eucl^Cka + Y (5) 

kcT q 

The symbol Ek is defined by: Ek = ~ whereas 

Ek is the electron band energy; p, denotes the chemical 
potential. We assume that electrons can propagate in 
the square lattice with the hopping integral t. In the 
considered case, energy band can be calculated using: 
Ek = —tj (k), where 7 (k) = 2 [cos {kx) + cos (ky)]. The 
function Wq determines the values of the phonon energies. 
The interaction operators are given by: 

= Y ^ (q) 4+q<TCkcT<?iq, (6) 

kqtr 


For T > Tc, the diagonal elements of Gk(ia;„) model 
the properties of the non-superconducting state. The 
condensate of Cooper pairs below the critical temper¬ 
ature is described by all elements of the matrix Green 
function. The Nambu spinors do not meet strictly 
anti-commutation relation [50] . Hence, the elements of 
Gk(*w„) should be calculated directly from the anti¬ 
commutation relation for the fermion operators. The 
equations for the scalar Green functions have a follow¬ 
ing form: 


(*w„) = (iw„) (za;„), 

( 10 ) 


gif'* (iuJn) = 5ok^ (*w„) mlf'’ (iuJn) ^ (*w„), (11) 


and 


470, E 

kk ql {7 


( 2 ) 
V , 

k.k 


' (q> 1) cLi<TCk<TcG 


k'-rn-q 


_„^k -o 


(7) 

where: = &)Lq -I- bq. The quantities (q) and 

(‘2'\ 

k k' denote the EPh and EEPh potentials, respec¬ 
tively. Next, we define the Nambu spinors: 


gY (*w„) = ^ (iun) mY (*w„) g^k^ (iuJn ), (12) 


5k^^ (ioJn) = 5ok^ (*‘^n)+5ok^ (*^") (*w„) ^ (iw„), 

(13) 

where the non-interacting Green functions are defined as: 


C-h)- (8) 

The matrix Green function Gk(*a;„) = 
takes the form: 


( ^k^lcjk^ Ck^|c_k4 


( 5k^^ (*w„) glf^ (iuJn) \ 

V ^ (*^") (iw„) j 




and 



(14) 

The scalar Dyson equations were achieved only for diag¬ 
onal Green function. The non-diagonal elements describ¬ 
ing the properties of the superconducting state meet the 
truncated Dyson equation due to the existence of the con¬ 
densate only in the case of non-zero pairing interaction. 
The self-energies rn[f^ (iujn)-m\^^ (itOn) can be written 
as: 


AA) 


(iuJn) = - (ui -I- V2 {n)f - EE << Ck-qt|Ck_q.^ 4 ’q\ 4 '-<i 


(15) 


m 


(B) 


(iuJn) = {vl + 2 viV 2 (n)) - EE Ck—q'l'|c_k-fq4, ^^iujrn 


(16) 


_ l!(i 

/33 2^ 


N 

L ki 


iOrn -I 


^ E << Cku|c-kit 


JV E '^-k2tl'3k2t 


E Ck-qt|C-k+q4 </>q|</>-q A-A- 









4 


= {vl + 2 viV 2 {n)) 4 51 XI « clk+q;|Ck-qt </'ql^-q »^{uJ^-u^) 

^ iJrn q 


/33 ^ 


iOrm ~^rT 


1 


-y 

N ^ 


« clkiTiy^ » 


1 


^ X << Ckstlc-ka; >>* 


k2 


X . 1 ^ _l_ r-t I I 


.t |J 


-'-k+qil'^k 




(17) 


,(£>) 


(ia;„) = -(?;i +?;2 (^^))XXX « c^k+qi|c-k+q; »iu^« ((>q|(('-q »i(c3„-c^„) ■ 


q 


( 18 ) 


New symbols occurring in the formulas (15 )-(|18|) were de¬ 
fined in the following way: (q) —t vi, i^k^j'^(li, I 2 ) —>■ 

and (n) = ^ X)ka (cLcktr^ The phonon propagator 
can be calculated in the non-interacting case: 

r‘^0 ^ 

<< ((>q|(()_q = —2 / dUfq (n) — 2 "——2^19) 

Jo LUi + iL 

— Pq (*Wz) ) 

where: w/ = (27r//3) I, Fq (fl) represents the phonon den¬ 
sity of states, and ojq is the maximum phonon frequency. 


Obtained results enable the derivation of the Dyson equa¬ 
tion for the matrix Green function: 


Gli(iUJn) — Gok(^^n)TGok(^^n).^k(^^n)Gok(^^n)i (20) 


where Gok(*w„) = (jWn'O ~ ^kTs)” denotes the propa¬ 
gator for the non-interacting electrons. The matrix of 
self-energy has the form: 


Mk(ia;„) = (ui -I- V 2 {n)f - 5 ^ rsGk-q (iujm) T^Pq (iw„ - iwm) 

CJ^nq 

(n)^ 1 

+ - « X ('^0 “ '^ 3 ) Gk-q {iiOni) (tq -I- Ta) Pq (lW„ - ZW^) 

CL»m.q 

(jlG 1 

+ - « X ('^0 + '^ 3 ) Gk-q {iiOni) (tq - T- 3 ) Pq (zW„ - ZW^) 

^ 0Jrn<l 


R3 X /V X t^ki 12 AT X [^*^2 


N ■ 

UJrm ki 


'm2Ul2 


^ X (^0 “ '^ 3 ) Gk^q (fw^ 3 ) (to + T 3 ) Pq (iUn - iUJ 7711 ^^ 777 , 3 ) 


q 

,,2 


fl3 X ZV X (*Wmi)] 21 TV X [^*"2 


N 

tOm-i ~G;mq ki 


'm 2 )\2l 


X ('^0 + ^ 3 ) Gk-q (to - To) Pq {iuJn - 777-1 ^^7712 ^^ 777,3 ) j 


( 21 ) 


where tq-to represent the Pauli matrix base and Next, the distribution of self-energy was postulated: 
the following dependence was used: = 

— ((A1|P))5 • The symbol * is the complex coupling. 

Mk{iuJn) = ZWn (1 - Zk (zw„)) To + Xk (*Wn)r3 ( 22 ) 
+ ipk{iuJ„)Ti +^^{i0Jn)T2. 
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On the other hand, the equation (|20[) can be written as: Substituting the expression (221 into (231 and calculating 


the overt form of the Green function, it was obtained: 


G-' {iWn) = G-^ (icc„) - Mk (lUJn) ■ (23) 


Gk {iUJn) = -i’k ^ [*Wn2’k Tq + (Xk (*W„) + ffk) T3 + (y3k (*Wn) Ti + (p^ (ioJn) T2] 


(24) 


where: I?k (*w„) = (w^^k (*w„))^ + (xk (i^n) +£k)^ + should be inserted into (21): 
<pI^ (iuJn) + ^k (*^ra)- In the next step, the formula (24) 


A^k (*W„) = - i {vi + V 2 {n)f ^ Wm^k-q (i^m) I^k-q (*^" “ "^0 

^ UJrn<^ 

- (^^1 + V 2 {n)f ^ ^ (Xk-q (*Wm) + Ek-q) -Dk(l;q (^^^m) -Pq {i^n “ iW^) '^3 
+ (f;! + 2uin2 (n)) ^ ^ (/Sk-q {i^m) Pq (*W„ - iOJm) Ti 

^ UJrn<\ 

X ~ ^ ^ [^k2 (^^m2) + ^V^k2 (^^m2)] ^k2 (^^"^ 2 ) 


X ^ ^ ^k—q ) -^k—q -fq ^^m2 

q 

+ (ui + 2win2 (n)) ^ ^ ^k-q (^‘^m) P'k-q Pq (*W„ - ZW„) T2 

^ CL)m,q 


+ 


4u| 


4 51 iiuJrm ) - iv^ki (*Wmi)] ^ki^ (^‘^rni ) 


N 

~^mq k 


X ^ 51 (*W„,J + *‘/ 5 k 2 (^Wms)] P’k^^ (*W^2) 


X ^ ( ^k—q (^^^3) -Pk—q (^^^3) -Pq (^^n ^^7713) ^ 2 - 


(25) 


Comparing with each other the formulas (22) and (25), complexity, the isotropic equations were taken into ac- 


we get the Eliashberg equations. Due to the numerical count, and we assume z ;2 —t ^2/2: 


Z„ = 1+ + y (n)) — ^P(n-TO)w^Z^—^DqXw), 


(26) 


Xn = -{yi + ^ (n)) ^ 51 ^ 51 -^q ^ + I'l™ - p] 


(27) 
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= {vf+ V 1 V 2 {n)) m) If (m) (28) 

^ m q 

+ )4 X/ P ip - mi - m 2 - m 3 + 1) f^^fm^VrUi 

^ mi~m3 

qi q2 qs 


where: = Z{iwn), Xn = xi'i-^n), fn = fiioJn), 

and Dq{m) = Z)q(iWm)- It should be noted that fn 
and fn take identical values. For this reason, the dis¬ 
cussion is limited to fn- The function Pin) is given 
by: P{n) = and v = PujoI^tt. The chemi¬ 

cal potential should be calculated from the equation: 
(n) = 1 —/_X (s) tanh ^ w is the half-width 

of the electron band, and the density of states has the 
form [51]: p{e) = 6 iln|^|, where &i = —0.04687t“^ and 
62 = 21.17796t. From the physical point of view, the 
function p (e) characterizes the van Hove singularity lo¬ 
cated at the Fermi level. In the cuprates it is related to 
the Cu 02 planes and is observed close to the Fermi en¬ 
ergy. For example, in YBa 2 Cu 408 the singularity exists 
about 19 meV below ep ad, the Y point in the Brillouin 
zone [52] . 


IV. THE ANALYSIS OF THE 
SUPERCONDUCTING AND THE ANOMALOUS 
NORMAL STATE 

The Eliashberg equations set was solved for 45 Mat- 
subara frequencies, which was related to taking into ac¬ 
count above 90 thousands of terms in the triple sum over 
Matsubara frequencies. The calculations used the mod¬ 
ified numerical methods originally designed to solve the 
classical Eliashberg equations with the electron-phonon 
pairing kernel [53] . The solutions of the Eliashberg equa¬ 
tions are symmetrical: Xm = X-m+i, and 

fm = f-m+ 1 - This property greatly simplifies the nu¬ 
merical calculations. The functions Zm, Xm, and fm are 
stable for the temperatures above Tg = 8.7 K. Results ob¬ 
tained for the superconducting state (T € (To,Tc}) are 
collected in the figures [^(A)-(C). In all cases it has been 
adopted vi = 150 meV and the increasing values of the 
potential V 2 from 0.75 eV to 5 eV were taken into account. 
The average number of the electrons was chosen in such 
a way that the order parameter would take its maximum 
value for T = Tg. It can be noted that obtained func¬ 
tions have characteristic lorenzian shape, formally simi¬ 
lar to the solutions obtained with the help of the classical 
Eliashberg equations [331 111 ] . The temperature depen¬ 
dence of the Eliashberg equations solutions can be most 
conveniently traced by analyzing the form of A„=i (T), 
Zn=i (T), and (T). In the Eliashberg formalism, 

with good approximation reproduces the physical 


value of the order parameter, Zn=i correctly determines 
the ratio of the electron effective mass to the electron 
band mass, and Xn=i well defines the direct renormal¬ 
ization of the electron band energy. The results are pre¬ 
sented in the figure [^ Typical curves characterizing the 
phase transition of the second type were obtained each 
time (Fig. [^(A)). This means that the complete descrip¬ 
tion of the disappearance of the superconducting state re¬ 
quires the precise consideration of the retardation effects 
(the explicit dependence of the Eliashberg equations so¬ 
lutions on Wn), the strong-coupling effects {Z), and the 
many-body effects (x). For case when they are omit¬ 
ted (as it was done for the toy model [311 HI]), appears 
the non-physical phase transition of the first order. On 
the other hand, the wave function renormalization factor 
and the energy shift function weakly depend on the tem¬ 
perature (Figs. (B) and (C)), which is the typically 
expected behavior of the considered quantities. It is also 
necessary to mention that growth of the EEPh interac¬ 
tion leads to the strong renormalization of the electron 
band mass and the band energy, which is caused by the 
high values of Zn=i and Xn=i- 

The figures |^(A)-(C) show plots of A„=i, Zn=i, and 
Xn=i in the dependence of the average number of elec¬ 
trons (doping). The most interesting results were ob¬ 
tained for the order parameter (Fig. [^(A)). Our model 
predicts that the dependence of A„=i on (n) for the low 
values of the EEPh potential {v 2 = 0.75 eV) looks very 
much like the dependence for the classical superconduct¬ 
ing state induced by the electron-phonon interaction [53] . 
Slight difference in the shape of the discussed function oc¬ 
curs only in the vicinity of (n) = 1 , where characteristic 
maximum can be observed. Situation drastically changes 
for the high values of the EEPh potential {v 2 = 5 eV). 
The order parameter vanishes for (n) < 1. It also plots 
the curve whose shape is observed for cuprates m even 
the characteristic plateau of the function A„=i ((n)) is 
evident for (n) G (0.8, 0.9) [55]. In the case of the wave 
function renormalization factor (Fig. (B)), its values 
increase together with an increasing V 2 - They are not, 
however, so heavily dependent on the doping as the order 
parameter. Let us note that the relatively weak depen¬ 
dence of Zn=i on (n) is confirmed by the experimental 
values of the electron effective mass obtained for the su¬ 
perconductors La 2 -xSrxCu 04 and YBa 2 Cu 30 y [56j. The 
figure (C) presents the plot of the energy shift func¬ 
tion. Results prove that the strongest renormalization of 
the electron band energy exists in the area of low (n) at 
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FIG. 1: (A) The order parameter, (B) the wave function renormalization factor, and (C) the energy shift function on imaginary 

axis. 



FIG. 2: (A) The maximum value of the order parameter, (B) the wave function renormalization factor, and (G) the energy 

shift function as a function of the temperature. 


the high values of V 2 - 

Interesting results concerning the superconducting 
state can be obtained for doping < n >= 0.95, in the case 
of which, the highest value of the critical temperature was 
determined {Tq = 39.4 K). By studying the dependence 
of the order parameter on the EEPh potential, we found 
that initially increasing V 2 causes the strong increase of 
A„=i. Nevertheless, above 1.25 eV, the order parameter 


again starts to decrease and disappear for z ;2 = 5 eV (Fig. 
|^(A)). From the physical point of view, this means that 
the highest value of the critical temperature is observed 
for the relatively weak EEPh interaction. It turns out 
that the cause of the decrease of the critical temperature 
in the range of the high values of V 2 is connected with 
significant increase of the electron effective mass and the 
strong renormalization of the electron band energy (Fig. 


























<n> 


<n> 


FIG. 3: (A) The maximum value of order parameter, (B) the wave function renormalization factor, and (C) the energy shift 

function as a function of doping. 


|4](B)and (C)). 

The dependence of the order parameter on the po¬ 
tential vi looks completely different. For any value of 
doping, for which the superconducting state is observed, 
increase in vi always causes the increase in (and 

therefore also increase in Tc). Sample results are pre¬ 
sented in the hgure (A). Additionally, the hgures 
(B) and (C) show the plots of the wave function renor¬ 
malization factor and the shift function. It can be seen 
that increasing potential vi does not result in the signif¬ 
icant increase in Zn=i, the decrease in absolute value of 
Xn=i is also observed. 

Referring to the results of above discussion, it can be 
assumed that the superconducting state with the maxi¬ 
mum critical temperature would rise in the case, in which 
EPh interaction will be strong enough, and EEPh inter¬ 
action will be relatively weak so that both the electron 
band mass and the electron band energy would not be 
renormalized in the significant manner. Taking into con¬ 
sideration the significant dependence of vi and V 2 on the 
distance between atoms building the crystal lattice, it is 
necessary to consider the possibility of carrying out the 
series of the high-pressure experiments analogously as 
it was done for lithium, calcium, and hydrogenated com¬ 
pounds [SMI]- 

The Eliashberg equations on the imaginary axis enable 
to determine the exact value of the critical temperature. 
They also allow (with the good approximation) the cal¬ 
culation of order parameter, electron effective mass, and 
energy shift function. However, imaginary axis formal¬ 
ism is useless in the analysis of the superconducting state 
in the non-static area (w 0). In order to determine the 
significance of dynamic effects, A„, and Xn should 
be analytically continued on real axis: A„ —> A(a;), 
Zn —> Z {uj), and Xn —t x(w). The easiest way to do 


it is using the formula [62] : 


X{uj) 


Pi + P2U} + ... + PrOJ^ ^ 

qi + q 2 UJ + ... + -I- w’’ ’ 


(29) 


where: X G {A,Z, y}. Parameters pj and qj are deter¬ 
mined in the numerical manner. In addition, we have 
adopted r = 6. 

The figures (A)-(C) present exemplary results ob¬ 
tained for the case: vi = 150 meV, V 2 = 5 eV, and 
(n) = 0.8, whereas the selected values of the temperature 
were taken into account. It can be seen that in the range 
of the low frequencies (w S (0, ^ 100) meV for A (w), and 
we (0, ~ 75) meV for Z (w) and y (w)), the solutions of 
the Eliashberg equations take the real values. From the 
physical point of view, this means the lack of the damp¬ 
ing effects (53]. In addition, in the considered frequency 
range, the real parts of solutions are poorly dependent on 
w, and it can be assumed that: Re [A (w)] ~ Re [A (0)], 
Re \Z (w)] ~ Re \Z (0)], and Re [y (w)] ~ Re [x (0)]. The 
situation radically changes for w ~ 110 meV (in the case 
of A (w)) and for w ^ 80 meV (in the case of Z (w) and 
xiuj)), where the appearance of the strong maximums 
and minimums can be observed. Let us note that above 
150 meV, the solutions of Eliashberg equations become 
extinguished. The increasing value of the temperature in¬ 
fluences very differently on the evolution of A (w), Z (w), 
and X (w)- It can be seen that the absolute values of the 
order parameter monotonically decrease with an increas¬ 
ing T. Interestingly, in the critical temperature even in¬ 
dicate non-zero values A (w) for w ~ 110 meV - these are 
the remains of the strong low-temperature maxima and 
minima. In the case of Z (u) and we practically 

do not observe the clear influence of the temperature on 
their low-frequency values, but their maxima and minima 
evolve differently - for the wave function renormalization 
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FIG. 4: (A) The maximum value of order parameter, (B) the wave function renormalization factor, and (C) the energy shift 

function as a function of V 2 . 
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FIG. 5: (A) The maximum value of order parameter, (B) the wave function renormalization factor, and (G) the band energy 

shift function as a function of Wi. 
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FIG. 6: The form of the order parameter (A), the wave function renormalization factor (B), and the energy shift function (C) 
on the real axis. 
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factor they noticeably decrease and for the energy shift 
function they increase. It should be emphasized that the 
observed values of Z (w) and x (w) are significantly higher 
than the values of the functions in question, if we take 
into account only the electron-phonon interaction. 

The physical values of the order parameter, the wave 
function renormalization factor, and the energy shift 
function should be calculated in the framework of the 
real axis Eliashberg formalism with the help of the equa¬ 
tion ([^ and the following formulas: 

Z{T) = Re [Z (w = 0, T)], (30) 

X(T) = Re[x(c^ = 0,r)]. 

The results obtained for the cases analogous as in the 
figure are collected in the figure As expected, the 
differences between data obtained in the framework of the 
imaginary axis and the real axis formalism are very slight, 
which proves that the analytical continuation procedure 
was properly conducted. 

As it was already mentioned, the real axis formalism al¬ 
lows to study the properties of the superconducting state 
outside the static area. It is possible then to give the an¬ 
swer to question: what can be the reason of the creation 
of the anomalous energy gap in the electron density of 
states, observed experimentally in cuprates [IMI]- The 
figure presents the form of N (w) determined for the 
case: vi = 150 meV, V 2 = 5 eV, and (n) = 0.8. It was 
found that in the range of the temperature from Tq to 
T* = 290.1 K, the energy gap induces in the electron 
density of states at the Fermi level. In particular, for 
temperatures from Tq to Tq, the half-width of the en¬ 
ergy gap changes slightly, while the clear decrease in the 
maximum value of N (w) may be noticed. Interestingly, 
in contrast to the shape of the electron density of states, 
associated with classical electron-phonon interaction, the 
curve N (w) is clearly non-symmetrical with respect to 
the Fermi level. The observed asymmetry results from 
the strong electron band energy renormalization, which 
is caused by the high values of x (^^)- It should be noted 
that theoretically determined asymmetry of the electron 


density of states is very clearly visible in the shape of 
the experimental curves m- The energy gap wears off 
slowly above the critical temperature. The first non-zero 
value of the function N (oj) on the Fermi level can be ob¬ 
served only at the temperature T* equal to 243.7 K. For 
T > T*, in the course of the electron density of states still 
visible is the remnant of the energy gap, which gradually 
disappears with the increasing temperature. Finally, the 
pseudogap disappears at T*. The presented model pre¬ 
dicts that gap in the electron density of states, in the full 
range of the temperature (even below Tc), is induced by 
Z (uj) and x (w)- The overt contribution of N (cj) lodged 
by the superconducting state is negligibly small. This 
can be easily checked determining the shape of the elec¬ 
tron density of states on the basis of the formula 0 , in 
which it was adopted (p (cu) = 0. 

V. SUMMARY 

In the presented work we have partially reproduced 
the properties of the superconducting and the anoma¬ 
lous normal state in cuprates. Considerations have been 
based on the Hamiltonian modeling the electron-phonon 
and the electron-electron-phonon interaction. The ther¬ 
modynamic equations have been derived in the frame¬ 
work of the Eliashberg formalism. We have shown that 
the experimental dependence of the order parameter on 
the doping can be qualitatively reproduced for the suffi¬ 
ciently large value of the EEPh potential. Then, it was 
proven that the Eliashberg equations correctly character¬ 
ize the type of the normal-superconducting state phase 
transition. In our opinion, the most interesting result 
was obtained during explaining the origin of the energy 
gap in the electron density of states. We have shown 
that the energy gap is induced by the diagonal elements 
of the matrix self-energy (Z and x); which are associ¬ 
ated with the anomalous normal state above Tc- The 
explicit contribution to the course of the electron density 
of states derived from superconducting state proves to be 
negligibly small. 
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